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ABSTRACT 

One of the two breakout reactions from the hot CNO (HCNO) cycle is 15 0(a, 7) 19 Ne, which 
at low temperatures depends strongly on the resonance strength of the 4.033 MeV state in 19 Ne. 
An experimental upper limit has been placed on its strength, but the lower limit on the resonance 
strength and thereby the astrophysical reaction rate is unconstrained experimentally. However, 
this breakout reaction is crucial to the thermonuclear runaway which causes type I X-ray bursts 
on accreting neutron stars. In this paper we exploit astronomical observations in an attempt 
to constrain the relevant nuclear physics and deduce a lower limit on the reaction rate. Our 
sensitivity study implies that if the rate were sufficiently small, accreting material would burn 
stably without bursts. The existence of type I X-ray bursts and superbursts consequently suggests 
a lower limit on the 15 0(a,7) 19 Ne reaction rate at low temperatures. 

Subject headings: X-rays: bursts — nuclear reactions, nucleosynthesis, abundances — stars: neutron 



1. Introduction 

Type I X-ray bursts were discovered indepen- 
dently by Grindlay et al. (1976) and Belian et al. 
(1976) and occur in matter-transfering binary sys- 
tems in which a neutron star accretes hydrogen 
and helium from an unevolved companion star 
(Woosley & Taam 1976; Joss 1977, 1978; Taam 
1980; Ayasli & Joss 1982; Lewin et al. 1995; Bild- 
sten 1998). The transferred matter is heated to 
1 — 2 x 10 s K while it slowly descends into the 
neutron star atmosphere as freshly infalling mate- 
rial continuously piles on top of it. At these tem- 
peratures hydrogen burns stably into helium at a 
constant and temperature-independent rate given 
by the half lives of two nuclei participating in the 
/3-limited HCNO cycle (Wallace & Woosley 1981). 
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where the decay constants A = In 2/7^/2 are given 
by the half lives of 14 and 15 0, and Xhcno is 
the total mass fraction of these two nuclei. The 
triple-a reaction concurrently increases the mass 
fraction of HCNO material while depleting helium 
(Bildsten 1998). The mass fraction of HCNO ma- 
terial is important because it regulates the rate at 
which hydrogen burns to helium, determining the 
initial composition and the characteristics of the 
burst (Fujimoto et al. 1981; Fushiki & Lamb 1987; 
Cumming 2003; Narayan & Heyl 2003). 

The mass fraction of HCNO material can de- 
crease through processing into heavier isotopes 
via the 15 0(a, 7) 19 Ne and 18 Ne(a,p) 21 Na break- 
out reactions (Wallace & Woosley 1981; Wiescher 
et al. 1999). This can happen either during the 
burst, when peak temperatures exceed 10 9 K, 
or during the quiescent, non-bursting phase. At 
temperatures characteristic of the quiescent phase 
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(T < 2 x 10 8 K) only the 15 0(a,7) 19 Ne reaction 
allows significant leakage out of the HCNO cycle 
(Hahn et al. 1996). In this case the leak of HCNO 
material into heavier isotopes depends on the ratio 
between the a capture rate of 15 O and its decay 

N A < av > pX a /A a 
K = 7 ) K 2 ) 

A15Q 

where < av > is the thermally averaged 15 0(a,7) 
reaction rate per particle pair, Na is the Avogadro 
constant, p is the stellar mass density X a is the 
4 He mass fraction and A a its atomic mass. 

What happens therefore depends on the rel- 
atively unknown 15 0(a,7) 19 Ne reaction rate: If 
the rate is comparable to the Caughlan & Fowler 
(1988) rate, then the leak will be important in 
reducing the HCNO mass fraction of Eq. [1] and 
therefore in decreasing the rate at which hydro- 
gen burns into helium as matter sinks. With long 
burst recurrence times of thousands of seconds 
or even days (Strohmayer & Bildsten 2003), even 
small reaction rates can be significant, allowing 
hydrogen to survive to great depths. In addition, 
a rate of this magnitude would deplete 15 O above 
the ignition point during the burst runaway, en- 
suring that little or no 15 O survives for the next 
burst and requiring the accretion of a fresh layer 
of material. 

However, the a width of the dominant 4.033 
MeV state in 19 Ne is not well known, so the rate 
could be lower. If the rate is lower, there will 
not be any significant leak out of the HCNO cy- 
cle, whence hydrogen would burn to helium at a 
faster rate, depleting hydrogen from the ignition 
zone. In addition, a lower rate would prevent or 
weaken the breakout via 15 0(a, 7) 19 Ne ensuring 
less energy generation during the initial stages of 
the burst (Taam & Picklum 1979; Ayasli & Joss 
1982; Hanawa & Fujimoto 1984; Wiescher et al. 
1999). The survival of a significant fraction of 15 
in the upper layers would increase the hydrogen 
burning prior to subsequent bursts, thus altering 
the composition at the time of the runaway. 

In this paper we show that a sufficiently small 
15 0(a,7) rate can even prevent subsequent run- 
aways, causing hydrogen and helium to burn sta- 
bly. Stable burning creates copious amounts of 
12 C and other A < 20 nuclei which have implica- 
tions for superbursts. Astronomical observations 
may lead to better constraints on this rate and mo- 



tivate future measurements of the 15 0(a,7) 19 Ne 
reaction rate. Such measurements would help re- 
fine current X-ray burst and superburst models. 

2. The ls O(a,7) 19 Ne reaction rate 

For temperatures below 6 x 10 8 K, the 15 0(a, 7) 
reaction rate is dominated by the resonant contri- 
bution of the 3/2+ state at 4.033 MeV in 19 Ne, 
which is a-unbound by 504 keV. Appreciable res- 
onant contributions also come from the states at 
4.379, 4.600, and 4.712 MeV. At higher tempera- 
tures contributions from higher-lying states must 
be included. Non-resonant capture also plays a 
role, but is only really important below 10 8 K. 

The rate was first calculated by Wagoner 
(1969) and later revised by Wallace & Woosley 
(1981) and Langanke et al. (1986). In this paper 
we provide a new estimate of the rate based on 
recent experiments and calculations. We calculate 
three rates, a best rate and upper and lower limits. 

The S factor for non-resonant capture was cal- 
culated by Langanke et al. (1986) within the 
framework of the direct capture model of Rolfs 
(1973), adopting the known reduced a widths of 
the bound mirror states in 19 F. A constant S fac- 
tor of 23 MeV b was adopted by Langanke et al. 
(1986) for the reaction rate calculation. The reac- 
tions 15 N+a and 15 0+a are isospin mirrors lead- 
ing to analog states in 19 F and 19 Ne. Therefore it 
is justified to assume that the reduced a widths 
in 19 F and 19 Ne are similar, which is supported 
by a comparison by Davids et al. (2003). The 
non-resonant contribution was also calculated by 
Dufour & Descouvemont (2000) using the genera- 
tor coordinate method. Despite the fact that the 
authors of this work state that the El transition 
probabilities are overestimated and that the re- 
sulting S factor should be considered an upper 
limit on the non-resonant contribution, they find 
an S factor only 85% as large as the Langanke et al. 
(1986) result. For the non-resonant contribution, 
we adopt the Dufour & Descouvemont (2000) cal- 
culation for our best rate, 1/3 of it for our lower 
limit, and use the Langanke et al. (1986) calcula- 
tion for our upper limit. 

The contributions of resonances are propor- 
tional to r„r 7 /r tota ; with r Q( - 7 ) the 0(7) partial 
width and r tota ; the total width of the state; for 
the states under consideration here, no other chan- 
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nels are open, so T total = r a + T 7 . In addition, 
for the two lowest-lying states, T a <C T 7 , and 
their strengths are directly proportional to their 
a widths T a . These widths are not known for the 
resonances of interest and cannot be measured di- 
rectly without an intense, radioactive 15 O beam 
at low energy. 

Several indirect measurements of T a for these 
resonances have been published. Mao et al. (1996) 
analyzed the a reduced width of the analog state 
of the 4.033 level in 19 F and studied carefully the 
dependence of the result on the model parameters. 
These calculations led to an a width of 11.5 /xeV. 
An alternative approach is the measurement of the 
a-decay branching ratio r Q /r tota ;. Two recent ex- 
periments, Davids et al. (2003) and Rehm et al. 
(2003), found 90% confidence level upper limits 
of 4.3 xl0~ 4 and 6xl0~ 4 , respectively. To deduce 
the a width from this ratio, the total width must 
also be known. A lifetime measurement by David- 
son & Roush (1973) and a Coulomb excitation ex- 
periment by Hackman et al. (2000) set 2a limits on 
the total width, 6.6 mcV < r tota; < 445 meV. The 
lifetime of the mirror state in 19 F corresponds to a 
total width of 73 meV. We adopt the alpha width 
of Mao et al. (1996) for our recommended rate. 
The upper limit of 130 /zeV for the a width is set 
by the measured 90% confidence level upper lim- 
its on limits on T to tai and the a-decay branching 
ratio. Considering the uncertainties in the deter- 
mination of reduced a widths from states in the 
analog nucleus far below the Coulomb barrier, we 
consider an a width of 3% of our recommended 
value, 345 neV, as a conservative lower limit. 

For the 19 Ne state at 4.379 MeV, Langanke 
et al. (1986) assumed a reduced a width of 0.02, 
leading to an a width of 700 [ieV. Three exper- 
iments report a-decay branching ratios for this 
state, 0.044 ± 0.032 (Magnus et al. 1990), 0.016 
± 0.05 (Rehm et al. 2003), and a 90% confidence 
level upper limit of 3.9xl0~ 3 (Davids et al. 2003). 
The first of these was based on poor statistics 
and represents an upper limit within 2a, while 
the energy resolution of the second measurement 
was not sufficient to resolve the 4.379 MeV level 
from nearby states, limiting the reliability of the 
extracted a-decay branching ratio. Therefore we 
adopt the upper limit of Davids et al. (2003). The 
lifetime experiment of Davidson & Roush (1973) 
set a la lower limit on the total width of 5.5 mcV. 



A shell model calculation using the USD interac- 
tion cited in Davids et al. (2003) found a radia- 
tive width of 458 mcV. Wilmes et al. (2002) have 
measured the a width of the 19 F analog state, 
from which Davids et al. (2003) have deduced a 
reduced a width of 0.0078t£ool 8 - Adopting the 
shell model calculation for the radiative width and 
the reduced a width from the analog state leads 
to our recommended a width of 190 /xeV. To cal- 
culate a lower limit on the a width we use the la 
lower limits on the total width and the reduced a 
width and find 1 /zeV. Our upper limit on the a 
width is based on the shell model calculation for 
the radiative width and the experimental upper 
limit on the a-decay branching ratio; using these 
values we calculate an upper limit of 2.3 mcV. 

Above 1 GK, the states at 4.600, 4.712, and 
5.092 MeV make appreciable contributions to the 
reaction rate. The parameters of these resonances 
were taken from Davids et al. (2003). Using these 
results we calculated three different rates of the 
15 0(a,7) 19 Ne reaction: a recommended rate, an 
upper limit and a lower limit. The rates are shown 
in Fig. 1. Model 1 represents the upper limit for 
the reaction rate, model 2 represents the recom- 
mended reaction rate, and model 3 represents a 
conservative lower limit. 




T[K] 



Fig. 1. — The three 15 0(a,7) 19 Ne reaction rates 
adopted in this work compared to the rate of 
Caughlan & Fowler (1988). 
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3. Computational model calculations 

We calculate three X-ray burst models (Fig. [3] ) 
that differ only in their different 15 O (a, 7) 19 Ne re- 
action rates. These models are calculated using 
a new general relativistic type I X-ray burst sim- 
ulation code that couples the general relativistic 
hydrodynamics code, AGILE (Liebendorfer et al. 
2002), with the generic nuclear reaction network 
of Hix & Thielemann (1999) using the operator- 
split method of Henyey et al. (1959). Specifi- 
cally, the calculation of a single time-step splits 
into four steps: 1) The nuclear reaction network 
is evolved in all zones assuming fixed density and 
temperature. 2) The change in nuclear binding 
energy enters AGILE as a source term and the hy- 
drodynamics is evolved. This step includes the 
self-consistent accretion of the rest mass. This 
step also calculates the turbulent velocity that 
determines the convective mixing of the compo- 
sition 3) The convective mixing of the composi- 
tion is evolved by solving a simple diffusion equa- 
tion based on the turbulent velocity and a mixing 
lenght 1 , A = 1cm (see Weaver et al. (1978)). 4) 
Finally, the composition is updated using the ad- 
vection that was calculated in step 1 using the 
second-order scheme described below. 

The code uses an arbitrarily relativistic and ar- 
bitrarily degenerate equation of state and includes 
radiative, conductive, and convective heat trans- 
port as described in Thorne (1977). We calcu- 
late the radiative opacities due to Thomson scat- 
tering and free-free absorption using the analytic 
formulations of Schatz et al. (1999). We use the 
same conductivity formulations for electron scat- 
tering on electrons, ions, phonons, and impurities 
as Brown (2000). 

AGILE solves the general relativistic equations 
in a spherically symmetric geometry on a comov- 
ing grid. Fig. [2] shows how the computational do- 
main is discretized into 103 log-ratioed grid zones 
with a column density 2 ranging from y = 1.2 x 10 6 

1 Since convection is extremely efficient, the exact value of 
A is unimportant. We note that setting it equal to the 
pressure scale height, which is typically done in stellar evo- 
lution calculations is unphysical, since the pressure scale 
height of the neutron star atmosphere is much larger than 
the convective zone itself. 

2 The relativistic column density is mass of a column above 
an area: y = f£_ r p^f where V = S /1-2GM/Rc 2 , so 
P ~ gy, where R is the neutron star radius, M is the 



g cm -2 to y = 3.5 x 10 9 g cm~ 2 . In this paper, 
AGILE does not employ its adaptive grid capability. 
The composition is advected using a second-order 
total variation diminishing upwind scheme with a 
van Leer limiter. 

The computational domain is bounded by a 
realistic core boundary interface (Brown 2003, 
2004) and a relativistically corrected grey atmo- 
sphere (see Thorne (1977) and Weiss et al. (2004)). 
Due to numerical constraints, previous work (e.g. 
Taam (1993); Taam et al. (1996); Fisker (2004); 
Fisker et al. (2004, 2005a,b); Woosley et al. (2004)) 
has set the outer boundary at P ~ 10 20 crgcm~ 3 
while assuming no nuclear burning and an opacity 
dominated by Thomson scattering at lower pres- 
sures. While the assumption of neglible nuclear 
burning is correct, the assumption of a constant 
opacity is off by up to 15%. The reason is that 
free-free absorption still contributes to the opac- 
ity at this pressure. Since the grey atmosphere 
model relates the opacity and the luminosity to the 
temperature and the temperature gradient, an in- 
accurate opacity at the boundary affects the tem- 
perature profile through the entire model. This 
approximation is perhaps acceptable at low ac- 
cretion rates or very high accretion rates since it 
does not change the overall bursting behavior of 
the system, but it is very critical at the transition 
rate between stable burning and unstable burning 
as stability is determined by a very fine balance 
between the cooling rate of the radiation trans- 
port and the heating rate of the nuclear reactions. 
The solution is to use a 4th order Runge-Kutta 
method and numerically integrate the hydrostatic 
heat and pressure equations (see Thorne (1977) 
and Weiss et al. (2004)) from the model boundary 
out to P = 10 18 erg cm -3 where the opacity finally 
becomes constant. 

Although heavier ashes with A > 100 predomi- 
nate, iff the burst peak temperature and the hy- 
drogen concentration at ignition are high (Schatz 
et al. 2001), Woosley et al. (2004) showed that a 
self-consistent X-ray burst simulation have compa- 
rably smaller peak temperatures and a smaller hy- 
drogen concentration at the ignition point. There- 
fore most of the ashes concentrate around A « 64 
with only a very small fraction of the nuclei react- 

ncutron star mass, p is the density, P is the pressure, and 
g = GM/YR? is the surface acceleration of gravity. 
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Fig. 2. — The rest mass in a zone as a function of 
column density. The computational domain covers 
about 7 pressure scale heights and comprises 103 
zones which are log-ratioed in rest mass. 

ing beyond this mass range. This means that we 
can use a smaller reaction network without chang- 
ing the underlying burst physics of the opacity, 
the equation of state, or the nuclear energy release 
rate. 

The nuclear reaction network used in this work 
includes 298 isotopes (see table [1]). All the con- 
necting particle reactions are taken from REA- 
CLIB (see Sakharuk et al. (2006)). This includes 
all isotopes between the valley of stability and the 
proton dripline up to 64 Ge. Here isotopes with j3 + - 
half lives > 1 day are considered "stable" on the 
timescale of the burst intervals, so their daughters 
are not included. Above the 64 Ge waiting point, 
only isotopes between the proton drip line and 
half lives less than 1 minute are included. This 
is because protons only capture on these high- 
Z isotopes during the burst's peak temperature 
which is only sustained for a few seconds. Weak 
rates up to Z = 32 are taken from Fuller et al. 
(1980, 1982a, b) and Langanke & Martmez-Pinedo 
(2001). Since only a small fraction of material is 
processed above Z = 32, it is a reasonable approx- 
imation to ignore neutrino losses from heavier iso- 
topes (Schatz et al. 1999). These considerations 
significantly reduce the size of the network which 
decreases the simulation run-timc. 

Since this paper investigates the transition be- 
tween stable and unstable burning, we can not as- 



sume an artifical iron-atmosphere for the initial 
setup. The reason is that an iron-atmosphere pro- 
genitor allows more H/Hc to accumulate which 
results in a very energetic burst that could heat 
the atmosphere to such a degree that the ther- 
mal inertia (see Taam (1993)) of the initial burst 
would ensure that H/He would keep burning at 
the same rate it was accreted and thus resulting 
in stable burning. As iron-atmospheres never oc- 
cur on accreting neutron stars in reality, the ini- 
tial model was made by accreting and comput- 
ing several bursts until a limit cycle equilibrium 
was reached. At this point, the composition of 
the burst ashes was copied to the bottom of the 
model-envelope after which the limit cycle equilib- 
rium was reachieved as the model adjusted bound- 
ary conditions to the new temperature profile re- 
sulting from the changed thermal conductivity in 
the copied shells. This creates a model that is 
physically independent of artificial starting con- 
ditions and which is approximately self-consistent 
with the chosen reaction library. At the same time 
this avoids the need to compute the evolution long 
enough to replace the entire model envelope which 
is shown in Fig. 2. We selected a global (rela- 
tivistically local) accretion rate of M = 10 17 g/s 
(0.09MEdd), similar to the accretion rates used by 
Taam et al. (1996) and Woosley et al. (2004). The 
accreted composition is assumed to be solar (An- 
ders & Grevesse 1989). The underlying fiducial 
NS core has a gravitional mass of 1.4M Q and an 
areal radius of 11.2 km. 

Fig. [3] shows an initial burst for all models 
which can be ascribed to the fact that a limit cycle 
equilibrium (see Lewin et al. (1995); Woosley et al. 
(2004)) has not been reached. The first burst of 
any model is usually "special" , yet whereas models 
1 and 2 both repeat this bursting behavior, model 
3 turns stable after the second burst fizzles. As 
model 3 turns out to describe an oscillating form of 
stable helium burning in a thin shell (see Woosley 
et al. (2004)) with a period of around 4 hours, 
we performed a convergence study of this model 
to determine the required number of zones to re- 
solve this burning shell. Fig. [3] shows that the 
oscillation is resolved in models having more than 
100 zones in the covered range of column depths 
whereas 75 zones or less is insufficient. In the next 
section, we divide the models into two categories, 
bursting and non-bursting. 
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Table 1 : The table shows the list of isotopes which 
describes the rp-process. See the main text for 
details. This reaction network has been used in the 
following works (Fisker 2004; Fisker et al. 2004, 
2005a,b). 
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Fig. 3. — Luminosity (as seen from an observer at 
infinity) as a function of (observer) time for the 
various models. Model 3 (lower limit) has been 
calculated at two different resolutions which both 
become stable after the first burst. 



3.1. Behavior in bursting models 

Models 1 and 2 include comparably high rates 
and exhibit unstable burning/bursting behavior. 
During the quiescent phase the 15 0(a, 7) 19 Ne re- 
action establishes an outflow from the HCNO cy- 
cle, thereby reducing Xhcno and ensuring that 
hydrogen burns slowly and thus survives to a great 
depth. 

For model 2 the thermonuclear instability crite- 
rion for a degenerate atmosphere (see e.g. Rakavy 
& Shaviv (1967)) is violated at y = 7.9 xlO 7 g 
cm~ 2 . One hundred seconds before the burst peak 
temperature is reached, the ratio R of eq. [2] is 
<~ 1%; R exceeds ~ 10% about ten seconds prior 
to the burst which is just before the 14 0(a,p) 17 F- 
reaction becomes competitive and establishes the 
ap-process (Wallace & Woosley 1981; Fisker et al. 
2004). After escaping the HCNO cycle the reac- 
tion flow proceeds through the heavier isotopes as 
described in Woosley et al. (2004). During the 
runaway, the 15 0(a, 7) 19 Ne reaction will eliminate 
most of the 15 0, since 15 is only formed via the 
slow 14 0(/3+,T 1/2 = 76s) 14 N(p, 7 ) 15 0, thus re- 
ducing Xhcno- This means that the conversion 
of hydrogen to helium for the subsequent burst 
will be delayed, as explained in section 1. 

The final CNO abundance would also be influ- 
enced by the amount of accreted CNO material; 
small values, which might arise from spallation in 
the atmosphere (Bildsten et al. 1992), have been 
shown to create chaotic bursting behavior (Taam 
1993). 

3.2. Behavior in non-bursting model 

The low rate of model 3 (c.f. Fig. 1) stops the 
bursts and causes a slowly oscillating luminosity 
as seen in Fig. 3. This oscillating behavior will be 
explained now. 

A lower 15 0(a, 7) 19 Ne reaction rate restricts 
the outflow from the HCNO cycle, thereby causing 
a significant build up of the HCNO mass fraction 
at an earlier point (see Fig. 4). 

The large Xhcno mass fraction in turn con- 
verts hydrogen into helium at an increasing rate, 
fully depleting hydrogen at y ~ 8 x 10 7 g cm~ 2 . 
The associated energy generation raises the tem- 
perature of the envelope to T ~ 3 x 10 8 K and in- 
creases the rate of the 3-a reaction, transforming 
the helium into 12 C. This is immediately trans- 
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Fig. 4.— The mass fraction of H, 4 He, and the 
summed mass fractions of A~12, 14-21 nuclei 
for the stable steady state burning obtained by 
using the lower estimate of the 15 0(a, 7) 19 Ne re- 
action rate as a function of column density. This 
snapshot is taken at a time where the envelope is 
hot and the helium burning front is at y m 8 x 10 7 
g cm~ 2 , its uppermost level. 

formed into 14 O and 15 O as long as hydrogen is 
present. The higher temperature moves the he- 
lium burning front further out at lower densities, 
destroying hydrogen through the HCNO cycle up 
to column densities of y ~ 6 x 10 7 g cm~ 2 . 

After the hydrogen has been destroyed, energy 
generation decreases and the layer cools as it is 
advected down. Since there is no hydrogen, the 3- 
a reaction creates 12 C, which may capture another 
a particle to become 16 O. 

As the temperature of the envelope is now 
lower, the 3-a reaction increases Xhcno at a 
slower rate, whence hydrogen is burned slower. 
This allows hydrogen to survive to a depth of 
about y <~ 1 x 10 8 g cm~ 2 until the increasing 
temperatures finally speed up the 3-a reaction to 
make sufficient HCNO material to burn hydrogen 



off again. As a result the burning front of the 
hydrogen slowly moves back and forth between 
i/~6x 10 7 g cm~ 2 and i/~lx 10 8 g cm~ 2 . 

Even while burning in this stable mode (as de- 
fined by a one-dimensional model with a lower esti- 
mate of the 15 0(a,7) 19 Ne reaction rate) the neu- 
tron star may burst due to other reasons which 
are not considered in this model, e.g., localized 
burning or changes in the local accretion rate 
due to accretion instabilities. In this case a low 
15 0(a, 7) 19 Ne reaction rate will cause a restricted 
outflow of 15 during the burst, thus preserving 
HCNO material for the next burst. This causes 
the HCNO cycle to run comparably faster until 
the subsequent burst, which means that more hy- 
drogen is depleted by stable burning, thereby in- 
creasing the observable a-parameter of the X-ray 
burst source. 

4. Implications for superbursts 

Supcrbursts have been observed from a num- 
ber of sources (Wijnands 2001; Kuulkcrs 2002; 
Cornelisse et al. 2002) and have been theorized 
to be caused by carbon burning in the neutron 
star ocean (Cumming & Bildsten 2001). The 
narrow parameter space in which superbursts are 
possible is given in Cumming & Bildsten (2001), 
but current calculations (Schatz et al. 1999, 2001; 
Woosley et al. 2004; Fisker 2004) have not been 
able to generate the required amounts of 12 C; 
in addition, these calculations leave residual 4 He 
which will stably burn the 12 C into 16 as matter 
sinks deeper into the ocean. 

However, Fig. 4 shows that the ashes of sta- 
ble burning almost solely comprise 12 C along with 
other isotopes, mainly 16 0, 19 F, 18 F and 18 0. 
This means that stable burning of this kind need 
not be maintained at all times, but can work for 
short periods of a few months, which are sufficient 
to generate the required carbon for the next su- 
per burst. Cumming & Macbeth (2004) estimate a 
duration of weeks for the cooling time of a super- 
burst, during which the ocean is still hot enough to 
prevent X-ray bursts from occuring. This means 
that continuously accreting matter burns stably 
during this time, whence a significant fraction of 
the nuclei shown in Fig. 4 could be generated. 

One might also speculate that the local accre- 
tion rate is not the same over the entire star, and 
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that high local accretion rates might cause ma- 
terial to burn stably. It is therefore important 
to establish the value of the critical accretion rate 
where unstable burning turns stable (see Rembges 
(1999) and Fisker et al. (2003)) as well as obtain a 
good model of how the accretion flow distributes 
itself on the star. 

5. Conclusions 

We have demonstrated the significance of the 
15 0(a, 7) 19 Ne reaction to X-ray burst behavior 
and suggested a lower limit on the rate which could 
and should be pursued experimentally. 

Long intermittent periods of stable burn- 
ing, which may last for months, have been ob- 
served from burst sources; a lower value of the 
15 0(a,7) 19 Ne reaction rate will help stabilize sta- 
ble burning. The reaction rate determines the 
value of the critical accretion rate. 

In addition we have shown a way to generate 
the necessary carbon for the superburst theory to 
be valid. We have also shown that several different 
kinds of A < 20 nuclei are generated; applying this 
to superburst calculations requires proper knowl- 
edge of the fusion rates between these nuclei. 

We thank H. Schatz and S. Austin for suggest- 
ing an investigation of this rate and A.N. Os- 
trowski for discussions. JLF, JG, and MW 
were supported by NSF-PFC grant PHY02-16783 
through the Joint Institute for Nuclear Astro- 
physics 3 . BD acknowledges support from the Nat- 
ural Sciences and Engineering Research Council 
of Canada. 

REFERENCES 

Anders, E. & Grevesse, M. 1989, Geochim. Cos- 
mochim. Acta, 53, 197 

Ayasli, S. & Joss, P. C. 1982, Astrophys. J., 256, 
637 

Belian, R. D., Conner, J. P., & Evans, W. D. 1976, 
Astrophys. J. Lett., 206, L135 

Bildsten, L. 1998, in The Many Faces of Neutron 
Stars, ed. R. Buccheri, J. van Paradijs, & M. A. 
Alpar (Kluwer), 419 

3 see http://www.JINAweb.org 



Bildsten, L., Salpeter, E. E., & Wasserman, I. 
1992, Astrophys. J., 384, 143 

Brown, E. F. 2000, Astrophys. J., 531, 988 

— . 2003, private communication 

— . 2004, Astrophys. J. Lett., 614, 57 

Caughlan, G. R. & Fowler, W. A. 1988, At. Data 
Nucl. Data Tables, 40, 283 

Cornelisse, R., Kuulkcrs, E., J. J. M. In't Zand, 
Verbunt, F., & Heise, J. 2002, Astron. Astro- 
phys., 382, 174 

Cumming, A. 2003, Thermonuclear X-ray bursts: 
Theory vs. Observations, astro-ph/0309626 

Cumming, A. & Bildsten, L. 2001, Astrophys. J. 
Lett., 559, L127 

Cumming, A. & Macbeth, J. 2004, Astrophys. J. 
Lett., 603, L37 

Davids, B., van den Berg, A. M., Dendooven, 
P., Flcurot, F., Hunyadi, M., de Huu, M. A., 
Siemssen, R. H., Wilschut, H. W., Wortche, 
H. J., Hernanz, M., Jose, J., Rehm, K. E., Wu- 
osmaa, A. H. & Segel, R. H. 2003, Phys. Rev., 
C67, 065808 

Davidson, J. & Roush, M. 1973, Nucl. Phys., 
A213, 332 

Dufour, M. & Descouvcmont, P. 2000, Nucl. Phys., 
A 672, 153. 

Fisker, J. L. 2004, PhD thesis, Univ. Basel 

Fisker, J. L., Hix, W. R., Liebendorfcr, M., & 
Thielemann, F.-K. 2003, Nucl. Phys., A718, 
614 

Fisker, J. L., Thielemann, F.-K., & Wiescher, M. 
2004, Astrophys. J. Lett., 608, L61 

Fisker, J. L., Brown, E., Liebendorfcr, M., Schatz, 
H., & Thielemann, F.-K. 2005a, Nucl. Phys., 
A758, 447 

Fisker, J. L., Brown, E., Liebendorfcr, M., Thiele- 
mann, F.-K., & Wiescher, M. 2005b, Nucl. 
Phys., A752, 604 

Fujimoto, M. Y., Hanawa, T., & Miyaji, S. 1981, 
Astrophys. J., 246, 267 



8 



Fuller, G. M., Fowler, W. A., & Newman, M. J. 
1980, Astrophys. J. Suppl., 42, 447 

— . 1982a, Astrophys. J., 252, 715 

— . 1982b, Astrophys. J. Suppl., 48, 279 

Fushiki, I. & Lamb, D. Q. 1987, Astrophys. J., 
323, L55 

Grindlay, J. E., Gursky, H., Schnopper, H., Parsig- 
nault, D. R., Hcisc, J., Brinkman, A. C, & 
Schrijver, J. 1976, Astrophys. J. Lett., 205, 
L127. 

Hackman, G., Austin, S. M., Glasmacher, T., Au- 
mann, T., Brown, B. A., Ibbotson, R. W., 
Miller, K., Pritychenko, B., Riley, L. A., 
Roeder, B., and Spears, E. 2000, Phys. Rev. 
C, 61, 052801(R). 

Hahn, K. I., Garca, A., Adelberger, E. G., Mag- 
nus, P. V., Bachcr, A. D., Bateman, N., Berg, 
G. P., Blackmon, J. C, Champagne, A. E., 
Davis, B., Howard, A. J., Liu, J., Lund, B., 
Mao, Z. Q., Markoff, D. M., Parker, P. D., 
Smith, M. S., Stephenson, E. J., Swartz, K. B., 
Utku, S., Vogelaar, R. B., & Yildiz, K. 1996, 
Phys. Rev., C54, 1999 

Hanawa, T. & Fujimoto, M. Y. 1984, Publ. As- 
tron. Soc. Japan, 36, 199 

Henycy, L. G., Wilets, L., B6hm, K. H., Lelevier, 
R., & Levee, R. D. 1959, Astrophys. J., 129, 
628 

Hix, W. R. & Thielemann, F.-K. 1999, J. Comput. 
Appl. Math., 109, 321 

Joss, P. C. 1977, Nature, 270, 310 

— . 1978, Astrophys. J. Lett., 225, L123 

Kuulkers, E. 2002, Astron. Astrophys., 383, L5 

Langanke, K., Wiescher, M., Fowler, W. A., & 
Gorres, J. 1986, Astrophys. J., 301, 629 

Langanke, K. & Martmez-Pinedo, G. 2001, Nucl. 
Phys., A673, 481 

Lewin, W. H. G., van Paradijs, J., & Taam, R. E. 
1995, in X-Ray Binaries, ed. W. H. G. Lewin, 
J. van Paradijs, & E. van den Heuvel (Cam- 
bridge Univ. Press), 175 



Liebendorfer, M., Rosswog, S., & Thielemann, F.- 
K. 2002, Astrophys. J. Suppl., 141, 229 

Magnus, P., Smith, M., Howard, A., & Parker, P. 
1990, Nucl. Phys., A506, 332 

Mao, Z. Q., Fortune, H. T., & Lacaze, A. G. 1996, 
Phys. Rev., C53, 1197 

Narayan, R. & Heyl, J. S. 2003, Astrophys. J., 599, 
419 

Rakavy, G. & Shaviv, G. 1967, Ap&SS, 1, 347 

Rchm, K., Wuosmaa, A., Jiang, C, Caggiano, J., 
Greene, J., Heinz, A., Henderson, D., Jansscns, 
R., Moore, E. F., Mukherjee, G., Pardo, R.C., 
Pennington, T., Schiffcr, J., Sicmsscn, R., Paul, 
M., Jisonna, L., & Segel, R.E. 2003, Phys. Rev., 
C67, 065809 

Rcmbgcs, F. 1999, PhD thesis, University of Basel 

Rolfs, C. 1973, Nucl. Phys., A217, 29 

Sakharuk, A., Elliot, T., Fisker, J. L., Hemingray, 
S., Kruizenga, A., Rauscher, T., Schatz, H., 
Smith, K., Thielemann, F.-K., & Wiescher, 
M. 2006, in Capture Gamma-Ray Spectroscopy 
and Related Topics, ed. A. Woehr & A. Apra- 
hamian 

Schatz, H., Aprahamian, A., Barnard, V., 
Bildsten, L., Cumming, A., Ouellette, M., 
Rauscher, T., Thielemann, F.-K., & Wiescher, 
M. 2001, Phys. Rev. Lett., 86, 3471 

Schatz, H., Bildsten, L., Cumming, A., & Wi- 
escher, M. 1999, Astrophys. J., 524, 1014 

Strohmaycr, T. E. & Bildsten, L. 2003, New Views 
of Thermonuclear Bursts, astro-ph/0301544 

Taam, R. E. 1980, Astrophys. J., 241, 358 

— . 1982, Astrophys. J., 258, 761 

— . 1993, Astrophys. J., 413, 324 

Taam, R. E. & Picklum, R. E. 1979, Astrophys. 
J., 233, 327 

Taam, R. E., Woosley, S. E., & Lamb, D. Q. 1996, 
Astrophys. J., 459, 271 

— . 1977, Astrophys. J., 212, 825 



9 



Wagoner, R. V. 1969, Astrophys. J. Suppl., 18, 
247. 

Wallace, R. K. & Woosley, S. E. 1981, Astrophys. 
J. Suppl., 45, 389 

Weaver, T. A., Zimmerman, G. B., & Woosley, 
S. E. 1978, Astrophys. J., 225, 1021 

Weiss, A., Hillebrandt, W., Thomas, H.-C, & Rit- 
ter, H. 2004, Cox & Giuli's Principles of Stellar 
Structure (Cambridge, UK: Cambridge Scien- 
tific Publishers) 

Wiescher, M., Gorres, J., & Schatz, H. 1999, J. 
Phys G: Nucl. Phys., 25, R133. 

Wilmes, S., Wilmes, V., Staudt, C, Mohr, P., & 
Hammer, J. W. 2002, Phys. Rev. C 66, 65802. 

Wijnands, R. 2001, Astrophys. J., 554, L59 

Woosley, S. E., Heger, A., Cumming, A., Hoffman, 
R. D., Pruet, J., Rauscher, T., Fisker, J. L., 
Schatz, H., Brown, B. A., & Wiescher, M. 2004, 
Astrophys. J. Suppl., 151, 75 

Woosley, S. E. & Taam, R. E. 1976, Nature, 263, 
101 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 



10 



